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Abstract 



' Kovtun, Son and Starinets have conjectured that the viscosity to entropy density 

. ratio rj/s is always bounded from below by a universal multiple of h i.e., h/{ATTkB) for all 



forms of matter. Mysteriously, the proposed viscosity bound appears to be saturated 



Qh! in all computations done whenever a supergravity dual is available. We consider the 

• near horizon limit of a stack of M2-branes in the grand canonical ensemble at finite 



R-charge densities, corresponding to non-zero angular momentum in the bulk. The 
corresponding four-dimensional R-charged black hole in Anti-de Sitter space provides a 
holographic dual in which various transport coefficients can be calculated. We find that 
the shear viscosity increases as soon as a background R-charge density is turned on. We 
numerically compute the few first corrections to the shear viscosity to entropy density 
ratio r]/s and surprisingly discover that up to fourth order all corrections originating 
from a non-zero chemical potential vanish, leaving the bound saturated. This is a 
sharp signal in favor of the saturation of the viscosity bound for event horizons even in 
the presence of some finite background field strength. We discuss implications of this 
observation for the conjectured bound. 
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Introduction 



In recent years, string theory/gauge theory dualities have been proven to be immensely in- 
valuable tools in understanding strong dynamics of certain quantum field theories. For any 
strong-weak type duality, computations that are possible on both sides of the duality are 
scarce. That is where taking various limits of both sides of the duality comes in handy, 
and may provide evidence for this in the past include the PP-wave limit of the AdS/CFT 
correspondence [Tj among others. The hydrodynamic limit could be one such simplifying 
yet non-trivial limit where even highly quantum mechanical systems behave simply and uni- 
versally. Hydrodynamics appears to be relevant to achieve a better understanding of the 
quark-gluon plasma (QGP) state in heavy ion collisions at RHIC. Studying near equilib- 
rium phenomena in a hot, strongly coupled QCD plasma is never easy, even on the lattice. 
Extracting retarded Green's functions from Euclidean lattice computations requires long 
Minkowski time separations, which calls for a large number of lattice points [2^ . Therefore 
using a dual description for QCD like theories in order to extract transport coefficients may 
find even practical use in studying near-equilibrium QCD physics. 

Sometime ago, Policastro, Son and Starinets proposed a prescription for calculating 
Minkowskian field theory Green's functions using the supergravity dual JT]. Since then 
an extensive study has been done where, using this prescription, various transport coef- 
ficients were calculated from the gravity side corresponding to the D-brane world-volume 
theories [3] as well as M2 and M5-brane theories in M-theory |4. In 5j, it was conjectured 
that the ratio of shear viscosity to entropy density is bounded from below 

- > -4— = 6.08 X lQ~^^K.s, (1) 
s ~ Anks 

where t] is the shear viscosity and s is the entropy density. The bound is well satisfied 
for weakly coupled systems, which could be understood intuitively by noticing that in a 
weakly coupled plasma, mean free flight time for the constituents is long. The puzzling 
feature shared by every transport coefficient calculation performed so far ^Bj is that for all 
known holographic duals to various supersymmetric gauge theories at finite temperature, 
the proposed bound appears to be saturated This suggests that, at the infinite coupling 
limit where the supergravity description is adequate, there exists some sort of universality 
in the hydrodynamic description of all of these field theories (their dual being some black 
hole in anti-de Sitter space) 6]. It was also emphasized |E] that the universal nature of the 
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ratio is connected to the universality of the black hole absorption cross section for low energy 
graviton scattering, at least for the cases where the AdS near horizon region has a flat space 
completion. 

Just on dimensional grounds, the bound itself appears to be in harmony with the observa- 
tion that ?7/s is a product of the energy per effective degree of freedom in the field theory at a 
t'Hooft coupling Qym-^ ^ 1 ^^"^ ^ time scale associated with the mean free flight time of the 
quasi-particle excitation [6^. According to the Heisenberg uncertainty principle, this product 
must be bounded from below by a multiple of h in order for the notion of quasi-particle to 
make sense. In jjj , the universal nature of the ratio was further established and extended to 
a large class of supergravity backgrounds where the dual possesses a translationally invariant 
horizon. The entire class saturates the bound. 

There is no known example where, at the infinite coupling limit, the bound holds but isn't 
saturated. One major motivation behind this work was to investigate whether there exist 
example (s) in which, at the strict infinite coupling limit, the viscosity bound is satisfied but 
not saturated. A particular setup which could potentially avoid the no-go theorems discussed 
in [2j , involves supersymmetric gauge theories living on the world- volume of type II D-branes 
or membrane and five-brane theories in M-theory at finite global charge densities. From a 
lower dimensional prospective this corresponds to black holes with some finite gauge field 
strength turned on at the horizon. It is natural to believe that the hydrodynamic properties 
of these horizons with a non-zero gauge field must be different from the " neutral horizons" . 

Here we work in the grand canonical description at finite chemical potential corresponding 
to finite R-charge i.e., 7^ where j° is the R-charge. Gravitationally, the R-charge arises 
from finite transverse rotation from a 10 or 11-dimensional point of view. 

We find that turning on a finite R-charge background increases the viscosity. In the 
conclusion section we speculate as to why this happens. To our surprise, we discover that 
the viscosity to entropy density ratio r]/s, remains the same as that for zero chemical potential 
case i.e., l/(47r) up to forth order in powers of Q/Th, where Q is the angular velocity and Th 
is the Hawking temperature. This provides clear evidence for the saturation of the bound 
even for the horizons with a finite gauge field. 

The paper is organized as follows. We begin with a lightning review of the hydrody- 
namic limit in systems with many degrees of freedom. After reviewing the prescription for 
computing Minkowskian retarded Green's functions in chapter 3. The relevant gravitational 
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background, i.e., R-charged Anti- de Sitter black hole is discussed in chapter 4. Gravita- 
tional perturbation theory of AdS4 black holes is the subject of discussioin in chapter 5. In 
chapter 6, We consider (T^yTpfj) correlation in Minkowskian signature using the prescription 
reviewed in chapter 3. Finally, we proceed to calculate the transport coefficient of interest 
here namely the shear viscosity rj and rj/sin chapter 7. 

2 Hydrodynamic Limit of Quantum Field Theories 

Linear response theory is the mathematical theory of the relaxation of small disturbances 
around equilibrium where the thermally averaged Minkowskian Green's functions (retarded, 
in order to account for causality) of the unperturbed system fully characterize the system's 
response to the external stimuli 

Gl...,,.M = ld'xe-^^-e{t){[Q,Ux),Q,,...m)p. (2) 

Here, G^_^ ^^...{q) denotes the retarded Green's function, Q^^,,Xx) is the operator correspond- 
ing to the conserved current which couples to the external world disturbance, /ijS are some 
spacetime indices and P is the inverse temperature. It is argued, and in simple cases explic- 
itly shown, that the slowly varying (both in space and time) behavior of the Minkowskian 
Green's functions of interacting field theories has a specific pole structure imposed by the 
"hydrodynamic equations" . These hydrodynamic limit conditions are usually satisfied when 
local thermal equilibrium is achieved. To create such circumstances, a fiuid must be in its 
high collision regime where interactions are important to the dynamics. Hydrodynamics 
is the study of small, long wavelength and low frequency fiuctuations of a medium in the 
vicinity of its equilibrium point^ In this limit all the fine structure of cutoff-scale physics 
gets wiped out, leaving only a few transport coefficients at low energies and long distances. 
It turns out that the relevant degrees of freedom to a hydrodynamical description are the 
charge densities of various global symmetries at the UV cutoff scale, along with the phase 
of the order parameters if any phase transition exists jH]. This is in accord with the fact 
that relaxation of any disturbance in conserved densities, in the deep IR limit, diverges. 
This makes charge densities the only relevant degrees of freedom in the hydrodynamic limit. 
Transport coefficients appearing in these sets of equations are not themselves part of the 
^Clearly, in a conformal theory at finite temperature, the temperature provides the required length scale 
for the system. 
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hydrodynamic description but, rather, are inputs. These coefficients could, in principle, be 
calculated from the slow varying part of the two point Green's functions like 0- It is exactly 
this type of computations which we will be concerned with in this paper. 

A fluid in equilibrium has a spatial energy- momentum tensor of the form Tij = p6ij. 
Slightly away from the equilibrium, extra stresses will be present as a result of viscous forces. 
The viscous part of the energy-momentum tensor of a fluid is proportional to the symmetric- 
traceless combinations of the momentum gradient of the fluid. It takes the following form 

T., = S.,p-v{^^+^)~iC-lv)S.,Vn, (3) 

where t] denotes the shear viscosity and microscopically is inversely proportional to the mean 
scattering rate. It is directly proportional to the mean free path for the effective degrees of 
freedom at a given value of the coupling. A strongly coupled fluid will have less viscosity 
as the mean free path is small and the energy in the perturbation taking the system away 
from its equilibrium, gets redistributed among degrees of freedom very quickly. This simply 
means that strongly coupled fluids are better approximations to the hypothetical notion of 
an ^^ideal fluid\ 

Also, C is the bulk viscosity, p is the local momentum density, and vr is the momentum 
flux. In a conformal theory, = 0. The energy-momentum tensor (jHl), along with the 
conservation equation for the energy-momentum tensor 

dtT"^' + dj{T' - p5'^) = 0, (4) 

atT°° + 9,T°' = 0, (5) 

form the complete set of the hydrodynamic equations specifying the system in its hydrody- 
namic limit. The natural frequency of the system specified by and (jSj), which we will 
be focusing on in this paper, is called the shear mode and is shown to possess the following 
dispersion relation 

u = -iVq^, (6) 

where V = ri/{e+p) is called diffusion constant, e is the internal energy density of the system 
and p is the pressure. Noticing that the dispersion relation (jHS)) appears as the pole structure 
of G^{q) in the hydrodynamic limit, enables us to read off the transport coefficient rj. 
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3 AdS/CFT in Minkowski Signature 



In order to explore the hydrodynamic limit, one has to find a way to calculate the two point 
functions in Minkowskian signature. The celebrated Anti-de Sitter/Conformal Field Theory 
duality is naturally formulated in Euclidean signature where the boundary configuration 
of the on-shell closed string background in the bulk of AdS acts as a source for generating 
all the Green's functions of the corresponding boundary operator. Formally, one can write 

{U[x,,)U[x,,)...U[x,J) oc 10 



where represents the closed string boundary value associated to the operator 0{x) 

in the boundary theory and Sd is the on-shell action. Calculating Minkowskian Green's 
functions is blocked due to the characteristics of the boundary value problem for hyper- 
bolic operators in Minkowskian signature AdS spacetime. The authors in pi] argued for a 
Minkowskian prescription for computing the Minkowskian thermal correlators in the bound- 
ary theory. According to this prescription 

Gniq) = -2T{u,q)\u=u,,^„,^,^, (8) 

where u is the radial coordinate in AdS to be defined later, q is the momentum on the 
boundary, and the imaginary part of is the Fourier component of the fiux associated with 
the corresponding AdS fiuctuation. The contribution at the horizon is discarded from (jH}. 
Moreover only incoming solutions to the AdS equations of motion are kept at the horizon. 
The prescription was further examined and confirmed to reproduce the desired and known 
results for the two examples worked out in JT] i.e., zero temperature A/" = 4, SYM theory 
in four dimensions and 2-dimensional finite temperature CFT dual to the BTZ black holes 
in AdSs. To see the detail of the definitions and justifications refer to jTT]. 

4 Gravitational Background 

4.1 D = A, M = 2 U(l)^ Extended Gauged Supergravities 

Extended gauged supergravities arise as Kaluza-Klein reductions of both D=10 and D=ll 
supergravities. Amongst these compactifications, an S'^ reduction of D=ll supergravity to 
N=8 SO (8) gauged supergravity in D=4 admits a consistent truncation where only gauge 
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fields in the Cartan subalgebra of tlie gauge group S0(8), i.e., a U{iy subgroup of the gauge 
group, survive. This truncation allows for 4-charge AdS black holes in four dimensions. Such 
a consistent truncation of D=4, N=8 S0(8) supergravity with f/(l)^ gauge group includes 
four commuting gauge fields, three dilatons, and three axions and is called minimal N=2, 
f/(l)^ supergravity. The Lagrangian for N=2 minimal supergravity is given by 

e"^£ = i?- i(90)2 + 8/(cosh0i + cosh02 + cosh03) (9) 



i=l 

where the 0j's are 3 dilatons, g is the inverse AdS4^ radius. The are introduced in |_12.^. It 
can be shown that any non-axionic solutions to this minimal supergravity can be uplifted to 
11-dimensional supergravity using the ansatz presented in |12j . 

4.2 4-Charge AdS4 Black Holes 

The near horizon geometry of the non-extremal M2-brane background which allows for up 
to four independent angular momenta (the dimension of the Cartan subalgebra of SO (8)) is 
given by [12] 



dsj, = K^'''[-{H,H2H^H,)-''^fde + {H,H2H^H,f'\r'dT^ + T^dy.dy)] (10) 
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where 



i=l 



f = -^ + Agh^H,H2H:,H,, (11) 
r 

= 1 + ^,^ = 1,2,3,4, 



X, = H-\H,H2H,H,y/\ 



gli sinh a ' 
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i=l 
1 



2 (2m sinh2a)"^/^f2, 
(2msinh^ a) 



r - 
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where the /i^'s are coordinates on the unit 3-sphere and h are four angular momentum 
parameters. The decoupled background (fTUI) gives rise to a duality between AdS^ x S'^ 
and the M2-brane theory at finite temperature and finite R-charge, such that in thermal 
equilibrium (j^) ^ 0, where j'^ represents the R-charge current in the boundary theory. 

The near horizon limit of the rotating M2-branes background under an S''' reduction 
gives rise to the four dimensional, 4-charge AdS black holes. As a result, it is more convenient 
to work within the framework of N=2 minimal supergravity described above while bearing 
in mind that any solution to this theory is an M-theory solution. Therefore we will focus 
only on constructing minimal supergravity perturbations. 

To simplify the setup without losing much of the generality, we choose to work with 
black holes with four equal charges which further simplifies the N=2 minimal supergravity 
to just Einstein gravity with a cosmological constant coupled to four Maxwell fields in four 
dimensions. In fact, setting four charges equal makes the dilaton sector decouple and partially 
simplifies the system without destroying its essential features. 

AdS4 black holes with four equal charges in this theory, with the horizon geometry being 
a space of constant curvature, are given by [12] 

dsl = -H-^fdt^ + H\f-^dr^ + r''dnlk), (12) 
r 

/isinh^/3 

n = L-\ ; , 

kr 

Ai = v^(l-if-^)coth/3,i= 1,2,3,4, 

where k = —1, 0, 1 refers to the curvature of the horizon geometry. 

The case k = needs special treatment and the result is the same background except 
the gauge field and the function H are now changed to 

H = i + (13) 
r 

The relation between the eleven and four dimensional metric is 

Ifg = 2/isinh2A, (14) 



■^Notc that in spacetimes with AdS4 asyniptotics, horizon topology is not restricted to just 2-spheres:the 
horizon manifold could be either of the 3 possible spaces of constant curvature. 
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smh/3j = glisinha, 



In what follows, we choose to work with the flat case k = 0. Upon uplifting to M-theory, 
this gives rise to the decoupling limit of the flat world- volume (parameterized by coordinates 
on ^2,0) rotating M2-branes with all four possible rotation parameters going. Note that the 
metric now is written as 

dsl = -H-'^fdt' + H\f-^dr^ + r\de^ + d<f)), (15) 

where 9 and are the dimensionless angular variables. ^ 

Let us introduce a more commonly used "u" coordinate for later use 



Rl 

r = — , 
u 

T^^ - JL 

^0 - 4^2' 



(16) 



where u is the new membrane radial coordinate. Note that functions / and H appearing in 
(fT^ are now written as 

/ = ^-^{H^-u'), (17) 
^ /i sinh^ (3u 

Here, we record a few quantities associated to the background (fTUI) for later usage. 

A dimensionless combination y = fx sinh^ P/Rq will make an appearance later on. The 
angular velocity Q corresponding to ()10j) is proportional to the chemical potential for the 
R-charge as viewed from the dual boundary theory. To calculate Q, we need to rewrite the 
M-theory embedding of our AdS^ black holes with four equal charges namely (jlUj) [T^ 

4 

dsl^ = -H~^fdt^ + H\f~^dr^ + r^dy.dy) + Y,idf^l + f^Kd<Pi + 9AldtY). (18) 



i=l 



Using definition of Q 



^. = _9t^^ (19) 



^The prescription for switching to the dimensionful coordinates x,y used in ^ is to notice that 2gx — 9 
and 2gy = (j). 
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where "i" labels each of the four independent angular velocities. One obtains 

"• = " = -»-4! = -*W- 

The event horizon is located where / = 0. The horizon radius can be expressed as a 
power series in sinh f3. We will keep terms only up to forth order in sinh f3 (or equivalently 



UH = Uo = l + ^y + 2y^ + 0{y'). (21) 



The Hawking temperature associated to this horizon is given by 
which leads to 



[22] 



Th = 3^(/i^^)V3(i_2 

TT 3 9 

where Tq is the Hawking temperature at = 0. The dimensionless ratio VL/Th can be 
expanded to the second order in y as well 

{^f = -^M^/^sinh/?(l + y+^y^ + 0(y^)), (24) 



-y{l + y + —y' + Oiy')) 



9 "9 

Using above and definitions (fT6|l . one can write y in terms of the ratio Q/Th up to forth 
order as follows 

y = A(^)'-2t?i(^)' + --- (25) 

5 Perturbing the R-Charged AdS4 Black Holes 

5.1 Review of the Reissner-Nordstrom Black Hole Perturbation Theory 

What follows is a lightning review of the perturbation theory of the 4-dimensional Einstein- 
Maxwell system with a cosmological constant. We closely follow [13 . Coordinates are 
labeled as {t,(j),r,6) = (0, 1,2,3). We follow the mostly minus signature convention for the 



spacetime metric. Perturbations are assumed to be generically non-stationary but axially 
symmetric. The most general non- stationary, axially symmetric perturbation of an arbitrary 
4-dimensional spacetime can be parameterized as follows 



dsl 



e^-dt" - e^^{d(P - udt - Qrdr - qedOf - e^^^dr^ - e^^'dO^ 



(26) 



It can be shown that the linearized perturbations fall into two distinct decoupled classes. 
One set, called "polar perturbations" , consists of 6Fq2, -fos, -^23, S^, Sfir, Sfie while the other set 
called "axial perturbations", includes Fqi, F12, Fi3,uj,qr,qe, where Fab denotes the Maxwell 
field strength. We use 6 in front of a fluctuation, whenever the corresponding fluctuation 
has a non zero background. We will not be considering polar perturbations here since the 
relevant perturbations to the viscosity computations fall into the axial perturbations class. 
The equations of motion governing perturbations are most easily written in the tetrad basis. 
The explicit form of the tetrad we use here is given by 



The hatted indices are fiat. Note that all the indices refer to the tetrad basis ()27|1 unless 
otherwise mentioned. 

5.2 The Axial Perturbation Equations of Motion 

The axial class of perturbation equations come from the following components of the Einstein 
and Maxwell's equations 

• (ab) = (12) and (13) components of the Einstein equations, 

• 1/ = (f) component of the Maxwell equations, 

• Bianchi identities written for (0, t, r) and (0, t, 6) permutations. 

The total number of equations sums up to five. There are two more equations which are 
redundant. In what follows, "," denotes ordinary derivative with respect to the corresponding 




(e^ 0,0,0), 

(-we^e^-g,e^-gee^). 



(0,0,e'^^0), 
(0,0,0,e'^«). 



(27) 
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coordinate. The explicit form of the axial equations of motion for the background (|T^ are 
written as follows 

{rf'/'Fo,),r+{H'rr'/'Fu),o = 0, (28) 
(r/^/^Foi),.+(i/VFi3),o = 0, (29) 

{H'rf-'/'Fo^),o+{rf/'F,2),r+{F,s),e = Hh^F,,Q,,. (30) 

where Qo^ = ~<lA,o and Q^b = QAjB ~<1b,a and A,B = 1, 2, 3. Taking derivatives with 
respect to r, 9 and t in equations (^Hj) . and pUj) yields 

(r/i/2Foi),,,,+(i/V/-V2Fi2),o„ = 0, (31) 

-fl3),0,6l — 0, (32) 

(i/V/-i/2Foi),o,o+(rr/'Fi2),.,o+(Fi3),e,o = i/VFo2Qo2,o • (33) 
From ()28|) we have 

[//7-2(r/V2FoO,.],.+(r/i/2Fi2,o),. = 0. (34) 
Utihzing dSH), (El, (El and dSl) we obtain 

[H-^f{rf/^Fo,),r],r+^{Foi)e,e-rH^r'^\Foi),o,o = -H^r^Fo2Qo2,2. (35) 

Now, let us turn to the Einstein equations. From (12) and (13) components of the Einstein 
equations, we obtain 

Ri2 = -V2'^-'^-'^«[(e3'^+'^-^'-^«g32),e-(e''^-'^-^'-+^«Qo2),o] (36) 

= —2Fr)TF'if\. 



Using (fT^ and (EHl) we get 



^2^3r/^32,e-i^VQo2,o] = -4FioFo2. (37) 
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The Ri3 component yields 



{r'fQ23),r = (i/V/-^)go3,o- (38) 



Equations (j?7j) and (jHHj) simplify to 



1 4 ^ 

jj^:iQ^e = —{^,r—qr,o),o+jp^f FiQFo2, (39) 

/ 

>Q,r = {^^,e —qe,o),o , 



where Q = r^f{qr,e —Qoyr )• Since dt and dg are killing directions of the unperturbed back- 
ground, a typical fluctuation will have the following form 

at,r,e) = ^{r)e'-'+'^', (40) 

where ^ denotes a typical fluctuation. ^ 
Eliminating uj from (j39p leads to 

9r[j^2Q^r] + j^,d^eQ-^d^Q = -^/"F,A,e. (41) 
Now let us return to equation Using (fHIIjl . one has 

[if-V(rf/'i^oi),.],.+^(Foi)e,.+(^V/fV-i/2_4^^2^^i^^ = -^Q,.. (42) 

Substituting the assumed form for the fluctuations pUj). one obtains the following pair of 
equations 

[i^-V(r/^/^FoO,.],.-^(FoO + (aWr^/2-4rFoV^/^)Foi = (43) 

It is also useful to work out the second order differential equation satisfied by the fluctuation 
UJ. In order to do so, the Einstein equation corresponding to i?oi needs to be written down 

Rm = -^e-2'/'-^'-'^«[(e3^-^-'^'-+^«g2o),r-+(e''^-"+^^-'^''g3o),e] (44) 
= 2F02-F12. 



^"g" appearing in 14()|l is dimensionless. Therefore to compare our results with one needs to send the 
dimensionless q to q/2g. 
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Using the AdS black hole background fields in ()12j) and simplifying the resulting expressions, 
one is led to 

I^{[H'r\qr,o-u;,r)],r+[H^r^f-\qe,o-u^,e)U) = -4^02^12- (45) 

Notice that utilizing the tetrad basis definitions in (jTfjl . the spacetime Fqi and Fq2 are given 
by 

J-oi = r/i/^Foi, (46) 

^02 = , 

where denotes the curved spacetime F. Using the convenient gauge where Qr = 0, (jl^ 
leads to 

Rewriting equation (^3)) in the gauge qr = gives 

-iH'r'u;,r),r+H'r'f-\qe,o,e-uJ,e,e) = -^^j^H'r\ (48) 

Simplifying the above equation utilizing definitions given in we will obtain 

{H^.)' -Ih^J -^^^^^{aqqe + q^.) = ^-^F,A, (49) 

where "prime" refers to d/du and A = Aq'^Rq. Using and (j^Hj) combined with well the 
gauge condition and (fT^ . one ends up with 

-^^{H^ - u')iqq,,^ = -iau;,,+^f'/^FoiFo2, (50) 

(r/V^FoO,. = taH'rr'/'F,2. 

Changing the coordinate system to "u" , equation (jTTj) is written as follows 

d .H^-v^ d 1 a^R^lP q\ . afismh(3I^ 
d^^^I^d^^''^ + A ^im - n3) - IP^^'' = A 

Denoting h = H^uj' and using (j49p and (j5(Jj) one ends up with a second order ODE for h 

which has the following form in "m" coordinates 

h - — ^ + -)h + - + --^ = 52 

^A''{m-U^Y ''^ Au\H^-U^yHml f ^ '^A(/74_^3)/^2^02^1'«' 

where, jFgi = —iaAi = —iaAfp using the fact that dt is a killing direction and remembering 
that we are considering axially symmetric perturbations. 
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6 Solving The Coupled System of ODEs 



6.1 Singularity Structure, Boundary Conditions 

The system of coupled differential equations (jSH) and (fS^ forms the fundamental set of 
equations to be solved. As is clear, these ODEs are singular at u = Uq where Uq is the 
horizon location. In order to isolate the singularity at u = uq, we substitute the following 
ansatz into the above ODEs 

= (uo-uyP{u), (53) 
h = {uo-uyF{u). 

The regularity condition, in addition to the incoming boundary condition for the fluctuations 
F{u) and P{u) at u = uq, will fix the values of 7 and u (as will be computed later) where 
Uq is the horizon radius. Substituting the ansatz (j^Hjl into and gives 

F{u)" +V{u)F{u)' + Q{u)F{u) = n{u)P{u) + S{u)P{u)' , 
P{u)" +U{u)P{u)' + V{u)P{u) = W{u)F{u), 

where 

Viu) = -( + ^ ^ + -)' (54) 

Uq — U 1 — U'^H * U 

v[v-\) V . (u^H-^)' 2, 2 2 (u^H-^)' 



^^"^ ^{uo-uy^uo-u^l- um-^ ^ v)^ u^^ ul- um-^' 



Siu) 



A{H^-u^) A\H*-u^y Au^m-u^ 
A^{H^-u^y' 



Av?{H^-u^yH^Rl f '^"'^ A{H^ - u^){uq - u)Rl 

2u2. 



A{m-u^)Rl^^' 



. 27 3u^-AH'^H H . 
Viu) ^ [^ + ^(^4.^ + 2^)+ ^^^^^0 



W(m) 



jj, sinh /3i?Q 
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As it was mentioned earlier, v and 7 can be computed by demanding regularity for functions 
F{u) and P{u) at m = Uq. v and 7 are thus given by the following expressions 

^ = v = ±^^(1 + + + 0(2/')). (55) 

Note that above expression can only be trusted to the second order in sinh/?. The minus 
sign corresponds to the the incoming boundary condition at the horizon and according to the 
prescription for calculating Minkowskian retarded Green's functions is the right boundary 
condition. 

6.2 Solving the System in Power Series, Domain of Convrgence 

In this subsection, we find the solution to the system ()54j) in a series expansion form around 
u = Uq. The aim will be to see if the radius of convergence of the series is large enough to 
include the point u = 0, where one is actually interested in calculating the pole structure of 
the Minkowskian Green's functions. As is obvious from the ODEs fl54j) . there exist 4 singular 
points: u = 0, u = uq, u = 00, and u = —R^/ (/i sinh^ Normally the radius of convergence 
of a series solution as viewed on the complex plane of u, extends all the way to the next 
neighboring singularity. For small values of f3, which is what we are considering here, the 
point u = —Rq/{^ sinh^ (3) will be well outside the convergence circle centered around u = uq 
and encompassing m = 0. So there appears to be no obstruction to continue the expansion 
to M = 0. Before presenting the solution, let us repackage our coefficients 

0( ^ - r 7(7-l) , 7 , {u^H-')' , 2, , 2 ,2 {u^H-')' 
^^""^ ~ \u,-uy^ Uo-u^l-u^H-'^ u^^ ul-um-^' ^^^^ 
g2 ^ 52^4 Axu^fiinh^p, 



' 2 ""^ sinh P o ™ sinh P 



nu) = -^[AQ'x——-^-8^x 



S{u) = 8x 



u sinh P 



HBf, • 



W{u) 



Rq sinh P 



where x = fi/R^, S = a/ (gy/x), Q = qj \fx. In the hydrodynamic limit, we will be interested 
only in expansions of the functions F and P at most to 3rd order in 5* and Q. To be precise; 
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we need to keep terms proportional to S and and nothing. This is because of the fact that 
diffusion phenomenon always involves two derivative with respect to the spatial dimensions, 
while there is only one derivative with respect to time. 

6.3 Numeric-Symbolic Solution 

In this section, we present our series solution to the coupled system of ODE's for the gravi- 
tational and gauge fluctuations. 

Let us digress for a moment and focus on how many integration constants one should 
expect in the solution. We have two second order ODE's, which means there are four 
integration constants. Two out of four are fixed by requiring the regularity condition for 
F{u) and P{u) at u — Uq. The remaining two integration constants get fixed by imposing 
boundary condition at the boundary of AdS, i.e., at u — 0. Starting from the series solution 
ansatz 

oo 

F{u) = ^/,(k-ko)\ (57) 

i=0 
oo 

i=0 

our plan will be to solve for /j = fi{S, Q, y) and pi — Pi{S, Q,y) wp to a desired order "N", 
as a function of the two remaining integration constants (which will turn out to be /o and 
po). y = (9/47r2)(fi/rff)2 - 2(81/167r^)(Q/rff)^ + e)((Q/7V)6) is the combination introduced 
in subsection (4.2). These series coefficients will be further expanded to the first order in S 
and second order in Q order which are the only relevant terms in the hydrodynamic limit 

MS,Q,y) = ^io(y) + ^ii{y)S + ^i2(y)Q'', (58) 
Pi{S,Q,y) = UM + Ila{y)S + Ui2{y)Q'. 

We further expand — ^kiiu) ^-nd Uki — Ilki{y) in powers of y 

^ki{y) = (t>kio + (t>kiiy + (t>ki2y'^ , (59) 

nfei(y) = TTfeio + TTfeiiy -Kkiiy"^- 

The interpretation of the indices is clear. In order to keep our notations simple, we have 
dropped the explicit dependence of fi and pi (and consequently all other expansion coeffi- 
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cients) on /o and po. Now all we are required to do will be to compute the coefficients (j)kii 
and TTkii. 

At this stage, we need to impose boundary conditions. To do this, we will have to use 
the perturbation equation (j^9|l . Remembering that 

h = H^uj' = {uo - uyPiu), (60) 

and taking the u ^ Q limit of (jiUI) give rise to 

{H^J)'\^^,--H'J\^^,-\{aqq's+q^uj') = 0, (61) 

where superscript "0" refers to the boundary values of the fluctuations qg and u at the 
boundary of the spacetime. Notice that the full solution for F{u) must go to zero at m = 
in order for F(u) to be a regular solution of PU]) " . Thus ()61|) can be rewritten as 

-(i7V)'U_o = F,o{S,Q,y)fo + Fn{S,Q,y)po = j{aqq'g + q'u;^), (62) 

where, we have substituted the series solution for F. We have also and taken into account 
and explicitly indicated the fact that, all the series coefficients are expressed as a function 
of S,Q,y as well as /o and po. Similarly, taking the m — limit of gives 

Piu)\^^o = Pio{S,Q,y)fo + Pn{S,Q,y)po = Al (63) 

where refers to the boundary value of Ai. 

Equations ()62|1 and (j63|l provide us with a system of two linear equations for two unknown 
/o and Po which fixes the integration constants /o and po in terms of the boundary values q^ 
and uj^. 

The diffusion constant denoted by V is the location of the Gtx,tx{(l,cr,y) po\e which is 
ultimately related to the shear viscosity through the relation 

V = (64) 

e + p ^ ^ 

The pole is given by the dispersion relation 

a = -Wq^. (65) 



"in our series solution F{u~{)) = ©(y"^), which is zero since we have only kept up to two powers of y at 
every step of our computations. 
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One can easily convince oneself that the location of the desired pole is the zero of the 
determinant of a 2 by 2 matrix F made out of Pio, Pu, -Fio and Fu 

Pio Pu 
Pio Fu 

The condition det{r) = gives rise to the following pole structure 

S = -tX(y)Q^ + Oiy'). (66) 

In order to compute V, we need to switch to the dimensionful quantities. Remembering that 
S = a/{gy/x), Q = q/{y/x), one can write as 

-^M2/)(77r^)', (67) 



(27rro)/3 ^"^\27rro)/3' 
where Tq denotes the Hawking temperature at zero chemical potential. ** Comparing to the 
dispersion relation (jH^I), one deduces 

- - H- 

In fact, A = X{y) is the quantity we calculate numerically as an expansion in y. 
T Calculating rj/s 

We proceed to calculate the shear viscosity using the relation (|64j) . Conformal symmetry 
implies e = 2p, leading to 

V = 3pl? = , (69) 

ZttIo 

where we have used V = 3\{y) / {2hTq) from the previous chapter. Note that in the grand 
canonical ensemble one has 

P = (70) 

where 'Em2 = E — TS — JQ is the Gibbs free energy for the spinning membrane. The Gibbs 
free energy for M2 brane with four angular momentum turned on, can be easily calculated 
jl4j . Expanding to fourth order in powers of Q/Th, one gets 

_ 9X4 J7 ^ 



"Note that in order to compare our results at leading order with T, we need to send q — s- q/{2g). 
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Similarly, the entropy density is given by 



1 dEM2 2V2vr2 9 X 40 ,2 27 , O ,4 



-N'^'T^l + - -T^A^Y + ••••]• (72) 



Assuming 



\{y) = \, + X,y + \,y^ + Oiy'), (73) 



the ratio rj/s becomes 



- = + + (74) 

s zvr 'in'^ 3 Aq Tf/ 

+^(-^-^^ + ^)A^ + ...) 

= pii+u^r+cA' +...). 

Itt Ih 1h 

As we mentioned in the previous chapter, A is the quantity we compute numerically. Thus, 
calculating A will provide us with the corrections to the ratio at zero chemical potential i.e., 
l/(4vr). 

Surprisingly, the following numerical analysis presented in the following, illustrates that 
the second and fourth order corrections to the rj/s ratio asymptotes to zero as we keep more 
and more terms in the Taylor expansion for F and A,^. This clearly signals a saturation of 
the bound even in the presence of a non zero chemical potential. 

The ratio at zero chemical potential has been calculated before to be 1/ (47r) j3] , therefore 
we expect to get A(?/ = 0) = 1/6 at the leading order. Of course, this only serve as a 
consistency check to assure us that the numerics have been done carefully. Below, we have 
plotted Aq = A(?/ = 0) versus A^, the number of terms in the Taylor expansions for F and A^ 
i.e., dilD 
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A at Q = vs. N 
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(2, the {VL/Th)'^ coefficient is the leading correction and tends to zero as it can be easily 
inferred from the figure 
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The {VL/ThY coefficient namely (^4 also runs to zero rather quickly as we increase N 



Next to Leading Order vs. N 




N 



For completeness, we have included a shear viscosity plot versus for N= 30. rjo is 

the shear viscosity for y = 



Viscosity vs. y 
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8 Conclusion and Outlook 

Our results satisfies the viscosity bound conjecture proposed by Policastro, Son and Starinets 
[5]. The leading finite chemical potential corrections to the viscosity itself turn out to be 
positive. The ratio rj/s remained unchanged up to fourth order in VL/Th signalling the sat- 
uration of the bound. An interesting question is why shear viscosity increases even though 
the system is literally at "infinite coupling". Another way of putting is to say, what could 
make an "infinitely coupled" system , "less infinitely coupled" (in order for the viscosity to 
increase)! One could speculate that the reason behind this enhancement in shear viscosity 
might lie in some screening effect for the color interaction mediator at finite chemical poten- 
tial. At non zero chemical potential i.e., when the number of various particle species carrying 
different R-charge is imbalanced, the "gluons" mass receives correction from the chemical 
potential which could result in screening . This screening of the color charge weakens the 
effective interactions in the plasma which ultimately leads to a bigger mean free path. When 
this paper was being written, I became aware of two other works [20], |2lj in preparation on 
AdSs with similar results. While for the membranes (which were studied in this paper), no 
gauge theory description exists, for the AdSs system there is a gauge theory living on the 
world-volume of D3-branes. Using conformal invariance of N= 4 SYM theory, only based 
on dimensional grounds, one could argue that the gluons mass receives corrections which 
are proportional to ^ where /i is the chemical potential A similar situation occurs in 
perturbative QCD. In this case, one could speculate with more confidence that the screening 
effect may be the true reason for an increase in the shear viscosity of a hot gauge theory 
plasma at finite chemical potential. 

It was argued and quantified by Karch [I5J that the conjectured viscosity bound is con- 
nected to Bousso's Generalized Covariant Entropy Bound (GCEB). Given such an interesting 
interrelation, one could reinterpret the viscosity bound as a non- gravitational and empirical 
window to the realm of quantum gravity. It turns out that the viscosity bound is exactly 
what matter is required to obey in order for gravity to modify the light-sheets (motion of 
the viscous fiuid results in a stress Tjj-generated curvature) to prevent the GCEB from a 
catastrophic violation. The current formulation of the GCEB suffers from a number of prob- 
lems including "the species problem" . The species problem is the simple statement that the 
^^/i has dimension mass. Note that, this effect is on top of the usual finite temperature corrections 
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entropy of a system of field(s) confined in a region of space can grow simply by increasing 
the number of particle species while keeping the total energy fixed. This would lead to a 
violation of the GCEB. An exciting question TB' would be to ask whether violating the 
GCEB through the species problem, for instance, would lead to a violation of the viscosity 
bound. To address this question one has to perform similar calculations as outlined in the 
present paper for supergravity duals to the gauge theories with large Nj |r7| (see also [TU] ) 
at "finite temperature". A zero temperature realization was considered in |19j, where the 
field theory corresponding to the localized D2-D6 intersection is an A/" = 4, c? = 3 super Yang 
Mills gauge theory coupled to Nq hypermultiplets in the fundamental of the gauge group, 
where Nq is the number of D6-branes. Nq/N2 is kept fixed while N2,Nq ^ 1. The full super- 
gravity background (i.e., D6-flavor branes including back-reaction) has been worked out in 
)18j . where the fact that uplifted D6-branes to M-theory has a Taub-NUT space component 
proves to be helpful. The non-extremal version of [IHj is not yet known and it seems like 
a daunting task to carry out. It would be interesting to find the background at least in 
the form of an expansion series in the vicinity of the horizon. Extensions of similar sets of 
computations would teach us a lot about whether there is a violation of the viscosity bound 
at large Nf. 

Needless to say, finding analytic solutions to the coupled differential equations here would 
be invaluable as it could reveal analytic structure of the viscosity as a function of VL/Th- 

There is no Gat in the statement of the bound. This simply indicates that the conjecture 
is an statement about quantum mechanical matter without any reference to gravity. So it 
is natural to expect the existence of a proof or counterexample for the bound which only 
involves weakly gravitating quantum physics. The relevance of gravity seems to be solely a 
consequence of the fact that in order to get down to the saturation limit of the bound, one is 
required to go to exceedingly high values of the coupling which, in the light of AdS/CFT, is 
mapped to the physics of highly gravitating objects, i.e., black holes in Anti de Sitter space. 
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